Collective motions in globally coupled tent maps with stochastic updating 
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We study a generalization of globally coupled maps, where the elements are updated with prob- 
ability p. When p is below a threshold p c , the collective motion vanishes and the system is the 
stationary state in the large size limit. We present the linear stability analysis. 
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I. INTRODUCTION 

The globally coupled maps (GCM) are introduced as a 
simple model capturing the essential features of nonlin- 
ear dynamical systems with many degrees of freedom [[lj . 
One of the most interesting phenomena seen in such sys- 
tems is the emergence of the collective motion ||, |], [|. 
The collective motion is characterized by a time depen- 
dence of the macroscopic variable in the large size limit 
§ | H| I H jL§ P, 0, |l| |l§. In this paper, 
we consider a variation of GCM to include asynchronous 
updating. 

The general form of GCM is written in the following 
way 



K 



N 



%t+i(i) = f{xt(i)) + 9{xt{i')) 



(1) 



where t represents discrete time steps, i specifies each 
element, K gives the coupling strength, and N is the 
system size. All elements are updated synchronously in 
the deterministic way through the mean field 



1 N 



(2) 



Since the collective motions have been studied anal yti- 
cally in globally coupled tent maps [|[ ||, 0, |[ ||, [l(], [lj] , 
|l2| , [ill , we specifically consider tent maps as follows: 



(3) 



f(x) = a I - - |x| 
g(x) ee f{x) -J. 



Here / is a constant determined from the average of f(x) 
over the the natural invariant measure of the map x i— » 
f(x), i.e., 



/ = / f{x)p*{x)dx , 



(4) 
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where p* (x) represents the the natural invariant density. 
By the above choice of g(x), h t = is a stationary so- 
lution for the large size limit (N — > oo) of ([[]), and the 
corresponding stationary distribution is proportional to 
p*{x). Equation (^) looks a little different from well- 
known form of GCM system 



x t+ i(i) = (1 - e)f{x t {£)) + ± £ f&W) , (5) 



which is obtained as a mean-field approximation for the 
coupled map lattice with diffusion coupling. In the case 
of the tent map system, however, the diffusion form (|^) 
is scaled into (0) Q. 

The collective motions in GCM ([!]) with (^) are clas- 
sified into two types according to the gradient a of the 
tent map. First, in the case of a < 1, the synchronized 
chaos is stable. In this case, the map x h- > f(x) has the 
stable fixed point 



(6) 



2a + 2 



Thus, p*(x) — S(x — a;*), i.e., g(x) — f(x) — /(x*). Since 
the gradient of f(x) is smaller than 1, the difference 
between any pair of elements diminishes. Thus all the 
elements behave identically after some initial transient. 
Here we concentrate on the long-term behavior and as- 
sume the system is one-cluster state. Then temporal evo- 
lution for xt and ht is obtained as follows 



Xt+l 

fh+i 



aK 



(1 + K)\x t \ 



2 2 + 2a 
a [x* — \x» + (1 + K)h t \] 



(7) 



When K is so large that a(l + K) > 1, the fixed point 
h t = is unstable. Then, the motion of the mean field is 
one-dimensional chaos, which obeys ([?]). 

Second, in the case of a > 1, all elements are fully 
desynchronized and behave as if they are mutually inde- 
pendent. Nevertheless, the fluctuation of the mean field 
dose not vanish in the large size limit 0. Thus, the sys- 
tem has a nontrivial collective motion I| . 

From a realistic viewpoint, however, the synchronous 
updating is not always plausible as models of real systems 
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FIG. 1: Bifurcation diagram for evolution of the mean field 
as function of the updating rate p in the case of a = 0.9, K = 
0.7,7V = 10 5 . For given values of p, the value of the mean field 
is plotted for 64 successive steps after 10 4 steps as transient. 



J!?], |ll| 19, [2C], [2ij| . In some cases, for example, neural 
networks, an independent choice of the times at which 
the state of a given element is updated should provide 
a better approximation. Abramson and Zanette have 
numerically found that, for globally coupled logistic map 
with completely stochastic updating, the fluctuation of 
the mean field can vanish in the large size limit |l8j |. In 
this paper, we study the stochastic updating model as 
follows 



xt+i(i) 



f(xt(i)) + Kh t with probability p 
x t (i) with probability (1 — p) 



.(8) 

At each time step, update the elements with probability 
p satisfying < p < 1. When the updating rate p is equal 
to 1, the model (j|) becomes the synchronous updating 
model (Q). On the other hand, when p decreases to 0, 
the model (|^) approaches the completely asynchronous 
updating model. Thus, the value of (1 — p) represents 
the strength of the asynchronousnism. The purpose of 
this paper is to investigate how the collective behavior 
changes when the updating probability p varies, mainly 
by the linear stability analysis of the stationary state in 
the large size limit. 



II. NUMERICAL RESULTS 

In this section, we present the numerical results for the 
model (|J). First, we examine the case of a < 1. When 
p = 1, the synchronized chaos is observed. On the other 
hand, when p < 1, some elements are updated and the 
others are fixed at each time step. Hence, even when a 
pair of elements have the same value at a moment, they 
can have different values at the next time. As a result, 
the synchronized chaos is broken. When p is near 1 , the 
synchronized state is blurred slightly. As p decreases, a 
sequence of bifurcations are seen (in Fig. Fll). It resembles 
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FIG. 2: Examples of the return maps of the mean field are 
shown (a = 1.625, K = 0.4). While (a), (b), and (c) are 
obtained by the direct numerical calculation of ^ with the 
size of N = 10 5 , (d), (e), and (f) are obtained by the numerical 
calculation of ( ]ll| ) with 2 18 grids, (a) and (d) are for p = 1, 
(b) and (e) are for p = 0.95, (c) and (f) are for p = 0.9. We 
plot 2000 points after 10 4 steps as transient. 



the period doubling cascade. However, there is a finite 
jump at this bifurcation in contrast of the usual period 
doubling bifurcation. This discontinuity is due to the fact 
that the map f(x) is piecewise linear. When p is smaller 
than a threshold value p c , all elements fall into the fixed 
point (Q) and the mean field h t becomes 0. Thus, the 
collective motion vanishes below the threshold p c . 

Second, we examine the case of a > 1, where the non- 
trivial collective motion is seen for p = 1. Figs. ||(a), 
||(b), and |2|(c) show the motions of the mean field for 
some values of p, with a = 1.625, K = 0.4, and N = 10 5 . 
The amplitude of motion of the mean field decreases as 
p decreases. 

It should be noted that the dynamics of each element 
is not deterministic due to the updating rule. Thus, the 
mean field value is blurred if the system size is finite. 
Even for p = 1, the finite size effect works on the motion 
of the mean field as an internal noise H . For that reason, 
it is useful to consider the large size limit. In the large 
size limit, the ensemble of the elements is characterized 
by its distribution. In the synchronous updating case, 
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the evolution of the distribution function p t (x) obeys the 
nonlinear Frobcnius-Pcrron equation 



Pt+i(x) = T[p t (x);h t ] 

_ p t (y+) + pt(y- 



(9) 



where T represents the Frobenius-Perron operator, and 
y + and y_ are the two preimage of x, i.e., x = f{y±) + 
Kht- Here the mean field is determined in the integral 
form: 



(10) 



ht=l g{x)p t (x)dx . 

In the stochastic updating case (p < 1) 
pt{x) is described as 



the evolution of 



pt+i(x) =p F\pt{x)\h t ] + (l-p)p t (x) 



Despite the stochastic updating, the distribution function 
evolves in the deterministic way. 

In order to calculate ([ll]) numerically, we approximate 
the distribution function pt{x) by dividing the relevant 
interval of x into m small intervals. The evolution of 
the distribution is described by the m x m transfer ma- 
trix which depends on time through the mean field h t . 
Here we construct the transfer matrix by applying the 



method by Binder and Campos p2| . Figs. |2j(d) 
and ||(f) show the the motion of the mean field calcu 



ated 



by this method with the parameter values corresponding 
to Figs. ||(a), |[t>), and ||(c), respectively. The direct 
calculations of (|) compare successfully with the results 
of except for the fluctuation due to the finite size 
effect. As is seen from Fig. |](f), when p is small, the col- 
lective motion vanishes like the case of a < 1. It should 
be noted that, in the stationary state, all elements arc 
still scattered and behave chaotically in contrast to the 
case of a < 1. The fluctuation of the mean field resides 
for finite size systems (Fig. ||(c)). 
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FIG. 3: Phase diagram of the collective motion for a = 0.9. 
The numerical results are obtained for TV = 10 4 . Circles and 
(11) crosses represent the cases of (h 2 ) 1 ^ 2 > 10~ 6 and (h 2 ) 1 ^ 2 < 
10~ 6 , respectively. Here the angle bracket denotes the average 
value of the enclosed quantity over 2000 time steps. The solid 
line represents the theoretical prediction (|14[). 



From ([lO|) and (|lj), we obtain the dynamics of the mean 
field obeys 



h t+1 = {-ap - apK + 1 - p) h t 



(13) 



Thus, the stationary state is stable when |— ap— apK +1 — 
p\ < 1. Consequently, the threshold value p c is estimated 
as 



Pc = 



a(l + K) + l 



(14) 



In the case of Fig. § p c = 0.790 • • •. The theoretical 
prediction agrees well with the numerical simulation (see 
also Fig. |). 

Second, we consider the case of a > 1. The stationary 
distribution function P*(x) is expanded into series of step 
functions as follows P, [13| 



III. LINEAR STABILITY ANALYSIS OF THE 
STATIONARY STATES 

The result of numerical simulation indicates that there 
exists a threshold value for updating rate p c . It is ob- 
served that the collective motion vanishes, and the sta- 
tionary state (with distribution pt(x) = p*{x)) is realized 
for p smaller than p c . In this section, we present the lin- 
ear stability analysis of the stationary state to estimate 
the value of p c . 

First, we consider the case of a < 1. This case is sim- 
pler, because all elements have the identical fixed value 
x*(> 0) in the stationary state. Considering a small per- 
turbation from a;», we assume that every element has a 
positive value. Since pt(x) for x < is 0, the evolution 
of pt{x) is rewritten as 



Pt+i(x) = - pt 
a 



a-2x + 2Kh t 
2a 



(l-p)p t (x) . (12) 



p*( x ) = c k o(x - x k ) 



k=l 

= / fe (o) 

= ^{(/^'(/(o))}- 1 



(15) 



where 6{x) is a step function: 1 for x > and for 
x < 0. Thus, X k and C k represent the locations and the 
heights of the steps in p*(x), respectively. We choose C\ 
to satisfy the normalization condition 



k=i 



CkX k = I dxp*(x) = 1 . 



(16) 



If there exists such k p that satisfies f kp (0) = and 
/ fe (0) ^ for Vfc < k p , the sum over k is taken from 
1 to k p . Otherwise, the sum is taken from 1 to oo. 

Before treating the case of stochastic updating, let us 
analyze how the stationary state is affected by adding 
external force with infinitesimal amplitude 6 for p = 1 
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and K = 0. Here we assume that the external force 
changes x$ for every element by the given amount. Thus, 
when the force Sq is applied at t = 0, the distribution 
function at t = 1 is expressed as 



pi(x;S ) = p*(x - S ) 



(17) 



In the limit of So — * 0, the response of the mean field 
after r steps is written as 



L T 5o 



(18) 



where L T is the linear coefficient for the response with 
delay of r steps From (|l5|), we calculate L T as fol- 
lows 



(19) 



fe=i 



When the temporal series of the external force is given 
as {St}, the mean held h t is obtained as 



h = L T 5t~ T 



(20) 



within the linear approximation. Introducing a cut-off d, 
the state at the time i can be described approximately 
by (i-dimensional vector 



(St-l,S t -2, S t -3: ■ ■ ■ S t -d) ■ 



(21) 



In order to obtain the stability condition accurately, we 
must take the limit of d — > oo. When the components of 
Vf is denoted as v\, we obtain 



ht = ^2 Li v\ 



(22) 



The next step is to consider the case of p = 1 and K ^ 0. 
In this case, the mean field coupling yields the feedback 
force. The influence of the feedback force is described as 



v 



t+i 



Kht 



(23) 



Taking (|22|) into account, the evolution of v\ is described 
as 



; '/ + 



i = H ' / " v 
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(24) 



where the matrix Jij is given by 



( KL X KL 2 KL 3 
1 
1 
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KL T \ 





(25) 
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FIG. 4: The stability condition. If all solutions /ii of (|26]) lie 
within the large (unit) circle, the stationary state is stable for 
p < 1 (p = 1). Thus, all /^i he on the left of the dotted line, 
there exists threshold value p c . 



The characteristic equation of the matrix (E3) is given as 



X d = K^UX 



d—i 



(26) 



which coincides with the results of Refs. |12|, |23[. The 
roots of (26) are denoted as fa. In the case of the syn- 
chronous updating (p = 1), if all pi lie within the unit 
circle in the complex plane, the stationary state is stable 
& 

Let us now return to the asynchronous updating case 
(p < 1). Here we define the vector v t to satisfy (p2|). 
Thus, v\ represents the contribution from the past per- 
turbations through i times of updating. Taking into ac- 
count that the elements are updated with probability p 
and fixed with probability 1 — p, we obtain 



v i+l = ^2\P J ij + O 1 - P) S ij] v l ■ 

The characteristic equation for (^) is given as 

d 



(27) 



A-l+p 
P 



A-l+p 
P 



(28) 



Comparing (P8[) with (J2q) , A in ( p6| ) is replaced with (A — 
1 +p)/p in (|28l), and thus the stability condition for (|27| ) 
becomes 



|1 — p + | < 1 (for all i) 



(29) 



This condition means all pi lie within the circle with 
center (1 — l/p, 0) and radius 1/p in the complex plane 
(Fig. |i|). In the limit p — > 0, the condition j29| ) becomes 
Re(pi) < 1. Consequently, if all /ij satisfy Re{pi) < 1, 
the system has the threshold p c . 

For example, we investigate a simple case a = (1 + 
\/5/2), where / 3 (0) = 0. Thus, fc p is 3 and X fc is periodic. 
In this case, the characteristic equation (p6|) can be solved 
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FIG. 5: Phase diagram of the collective motion for a = 
(1 + \/~E)/2. The numerical results are obtained by calculation 
of (Jill) with 2 18 grids. Circles and crosses represent the cases 



of (J?) 1/2 > 1(T 12 and {h 2 ) 1/2 < KT 12 , respectively, 
solid line represents the theoretical prediction 



The 



Figure [| shows the correspondence between the above es- 
timation for p c and the numerical results. It indicates the 
good agreement, except for the lower right corner. The 
estimation implies p c tends to 1 in the limit of K — > 0, 
and we think that the gap in the corner appeared because 
the amplitude of the collective motion for small K is very 
small (estimated at exp(— C/K) for the case with p = 1 



small (e 

0H) 



For such special values of a, where {Xk} falls on a pe- 
riodic orbit, we can solve the characteristic equation by 
the above method. In this case, the characteristic equa- 
tion has no solution p, which satisfies Re(fi) > 1 for ade- 
quately weak coupling. Thus there exists the threshold. 
For the general values of a, however, it is difficult to solve 
the characteristic equation and we have not obtained the 
stability condition explicitly at present. 



analytically in the limit of d — > oo. From ( |15[ ) and (19) 
the linear coefficient {L t } is given as 



£371+1 



Ll 

L 3 



V5 



10 /F 

5 - V5 



(30) 
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Let us assume that /i satisfies condition \a\ > 1. Then 
we rewrite the characteristic equation ( p6| ) in the limit 
d — > 00 as follows: 



1 =KY J UpT l 



K 



From this, [i is given as 



Li/i 1 + L 2 (J, 



(31) 



1-M" 



KL 1 -\± yj (KLi - l) 2 



(32) 



IV. SUMMARY AND REMARKS 



This study have explored the globally coupled tent 
maps with stochastic updating. We introduced the up- 
dating rate p and examined how the collective behavior 
changes as p varies. In the case of a < 1, the collective 
motion has a sequence of bifurcations, which is similar 
to the period doubling cascade. In the case of a > 1, 
the amplitude of the collective motion decreases as p de- 
creases, For the both cases, we observed the threshold 
p c , below which the collective motion vanishes. We esti- 
mated successfully the threshold p c by the linear stability 
analysis of the stationary state. 

For weak coupling, the threshold p c is likely to remain 
near 1 as is seen for the example a = (1 + v5/2). Thus, 
a tiny asynchronousnism may extinguish the collective 
motion. Therefore, when GCM is used as model of real 
systems, we keep in mind that even if the updating rule is 
almost synchronous, the effect of asynchronous updating 
should not be ignored. 



When (45 - 2lV5)/2 < K < (5 + 3\/5)/2, the part in 
the square root in ( |32| ) is negative and thus the solution 
( |32|) is c omplex conjugate pair. Then, we obtain = 
VKTT. Thus, \n\ > 1 holds for < K < (5 + 3y^)/2. 
In this case, the threshold p c is given by 



Pc 



KLi - 3 



K 



(33) 
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